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Abstract 



The Somos 4 sequences are a family of sequences satisfying a fourth order bilinear 
recurrence relation. In recent work, one of us has proved that the general term in 
such sequences can be expressed in terms of the Weierstrass sigma function for an 
associated elliptic curve. Here we derive the analogous family of sequences associated 
with an hyperelliptic curve of genus two defined by the affine model y 2 = 4x 5 + c^x 4 + 
. . . + c\X + Co- We show that the sequences associated with such curves satisfy bilinear 
recurrences of order 8. The proof requires an addition formula which involves the 
genus two Kleinian sigma function with its argument shifted by the Abelian image of 
the reduced divisor of a single point on the curve. The genus two recurrences are related 
to a Backhand transformation (BT) for an integrable Hamiltonian system, namely the 
discrete case (ii) Henon-Heiles system. 



1 Introduction 

In recent work |24j . one of us has considered fourth order quadratic recurrences of the form 



where a and (3 are constant parameters. Such recurrences arise in the theory of elliptic 
divisibility sequences [H3 EHl EI] and their generalizations, the Somos 4 sequences (HHl EH] • 
In that context, both the parameters a, j3 and the iterates r n are integers, or more generally 
take values in Q or a Galois extension, and in that case the sequences have applications in 
number theory, as they provide a potential source of large prime numbers P^ ll5j. Moreover 
the Somos 4 sequences, denned by a recurrence of the form (jl.lj) , provide a simple example 
of the Laurent phenomenon: taking the initial data tq,ti,T2,t^ and the parameters a, (3 as 
variables, all subsequent terms r n for n > 4 in the sequence are Laurent polynomials in these 
variables. Fomin and Zelevinsky have proved that this remarkable "Laurentness" property 
is shared by a variety of other recurrences in one and more dimensions, with applications in 
combinatorics and commutative algebra (see JB] and references). 
In [23] the following theorem was proved: 
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Theorem 1. The general solution of the quadratic recurrence relation M.l\) takes the form 



r n = AB"°-^^>, (1.2) 
a{K) n 

where k and Zq are non-zero complex numbers, the constants A and B are given by 

A= ^_ 8 = ^ (1 . 3) 

and a denotes the Weierstrass sigma function of an associated elliptic curve 

y 2 = 4x 3 - g 2 x - 3 . (1.4) 

The values k, z and the invariants g 2 , g% are precisely determined from the initial data 
To? Tii r 2? r 3 and the parameters a, j3. 

In the next section we summarize some facts about elliptic divisibility sequences, Somos 
4 sequences and the details of the above theorem. In particular we explain how the result of 
Theorem 1 is connected to the second order solvable mapping 

fn+l = -T^J- U + £-) (1.5) 

Jn—ljn \ Jn / 

which is a degenerate case of the type of mapping studied by Quispel, Roberts and Thompson 
|37| . (See also [3EJ for some recent work on the global behaviour of real- valued solutions of 
such mappings.) 

Remark. The case a = 0, which was excluded from the statement of the Theorem 
in |21], corresponds to k being a half period, so that jp(«) is a branch point of E, but 
then the formula for r n is has a rather trivial alternating form: r 2 k = t (t 2 / r o) k /3 fc( - fc_1 ^ 2 , 
T 2 k+i — r i( r 3/ r i) k P k( " k ~ 1 ^ 2 - The map (J1.5)) is the autonomous version of the discrete Painleve 
I equation (qdPI) 

fn + i = ( aq n + l-V (1.6) 



fn— if n \ fn 

which has a continuum limit to the first Painleve equation [H21 EH] • The qdPI map (jl.fjjl has 
tau-functions that yield a sequence of g-polynomials (2B] , and in the autonomous case q — 1 
this map reduces to (jl.5j) . Matsutani has constructed some particular solutions of (jl.5j) using 
elliptic functions, and has also considered certain higher order recurrences associated with 
genus two hyperelliptic functions [22] (see section 3). The above Theorem guarantees that 
only elliptic functions are necessary to specify the general solution of the second order map 

The result of Theorem 1 can also be understood via the addition formula 

a(z + k)g(z — k) , , , . . 

„wUo' = P{K) ~ p(z) ( ] 

for the forward and backward shifted Weierstrass sigma function in terms of the p func- 
tion (see e.g. jH])- I n section 3 this leads us to derive a higher order generalization of 
the recurrence relation (jl.lj) by considering a suitable addition formula for the Kleinian hy- 
perelliptic sigma function associated to a curve of genus two. The hyperelliptic Kleinian 
sigma functions are a natural extension of Weierstrass elliptic functions to the case of higher 
genus (see e.g. |21 IE] and references). The addition formula we consider is a special case 
of the generalized Frobenius-Stickelberger formula in OEH], which is the exact genus two 



analogue of (jl.7j) . The main result of our considerations is to derive an eighth order bilinear 
recurrence whose terms are given by an analogue of the formula (J1.2)) . As a corollary, we 
also derive the solution of a family of sixth order nonlinear difference equations in terms of 
Kleinian p functions. The fourth section explains how this recurrence is related to the BT 
(integrable discretization) of a Hamiltonian system with two degrees of freedom, namely the 
integrable case (ii) Henon-Heiles system; this BT first appeared in [UI22]> an< ^ was P u ^ m 
an algebro-geometric setting in j2Hj. The extension to higher genus is briefly discussed in 
our concluding section. 



2 Elliptic divisibility and Somos 4 sequences 

The sequence 

0,1,-1,-1,-1,2,1,-3,5,7,... (2.1) 
is an example of an elliptic divisibility sequence. It is obtained from the recurrence 

Tn+lTn-l + (r n ) 2 

Tn+2 = , (2.2) 

T n-2 

with initial data taken as 

Ti = 1, r 2 = r 3 = r 4 = -1. 

The sequence can be consistently extended backwards for negative n G Z, to give an an- 
tisymmetric sequence with r_ n = — r n . Remarkably, despite the division by r n _ 2 at each 
iteration of (|2.2|) . the subsequent terms of the sequence are all integers, and they satisfy the 
divisibility property 

T n \ T m whenever n\m. (2-3) 

More generally Morgan Ward |4"5tl4*6] introduced a family of such antisymmetric sequences 
defined by recurrences of the form 

r n+2 r n _2 = (r 2 ) 2 r n+ ir n _i - nr^Tn) 2 , (2.4) 

which are derived by considering sequences of rational points nP on an elliptic curve E over 
Q. To obtain integer sequences of this kind it is required that 

r = 0, n — 1, t 2 ,t 3 ,t a <eZ with r 2 |r 4 . 

Using the addition law on E and considering the multiples nP = (x n , y n ) of a single point 
P, Ward derived the bilinear recurrence (|2.4|) for r n , with the general term being written in 
terms of the sigma function associated with the curve E, as 

Tn = (2-5) 

(See also PHi]].) 

Using the addition formula (jl.7|) for the Weierstrass sigma function, it is a simple exercise 
to use the formula ()2.5|) in order show that the terms of the elliptic divisibility sequence satisfy 
the Hankel determinant relation 



i n+m ' n—m 



(2.6) 



for all m,neZ. Starting from the Hankel determinant formula, it is then easy to prove by 
induction that all r n are integers with the divisibility property (|2.3j) . 



If we consider the same recurrence (|2.2|) but instead take initial data tq = T\ = t 2 = T3 = 
1, then we find the sequence of integers 



1,1,1,1,2,3,7,23,59,314,..., (2.7) 

known as the Somos 4 sequence (see [2II1II2])- (In fact this Somos sequence is just obtained 
by selecting the odd index terms of the elliptic divisibility sequence (|2.1|) . up to an alternating 
sign.) More generally, following the terminology of [23112], we refer to any sequence defined 
by a bilinear recurrence of the form (jl.lj) as a Somos 4 sequence, while the particular sequence 
above is denoted Somos (4). It turns out that any such sequence is associated to a sequence 
of points Pq + nP on an associated elliptic curve E: this fact was proved by algebraic means 
in the thesis of Swart |43j . which refers to unpublished results established independently 
by both Nelson Stephens and Noam Elkies. In [21], one of us gave an alternative complex 
analytic proof, leading to the construction of the functional form (jl.2|) of the general term, 
as in Theorem 1 above. 

The approach taken in [21] was to regard equation (jl.l|) as the bilinear form of an 
integrable map, analogous to the bilinear equation satisfied by the tau function for a soliton 
equation and then solve the initial value problem for the bilinear equation with 

specified initial data r , T\, r 2 , r 3 and parameters a, (3. The quantity r n may be regarded as 
being the tau function for the second order nonlinear map (jl.5|) . to which it is related by the 
substitution 

„ T n+l T n-l 



The map f)l ,5|) has a first integral, given by 

J := J{f n -1, fn) = fn-lf n + « ( h — J + ^ = J (fn, fn+l)- 

\Jn—l In J Jn—ljn 

The algebraic formula for J itself implies that the pair (f n -i, fn) lies on an elliptic curve for 
all n. In fact the general solution of the recurrence can be written in terms of the Weierstrass 
p function for a curve in the canonical form (jl.4j) . as 

fn = - p{z + nK). 

The construction of this curve E and the points Pq, P E E solves the initial value problem 
for (jl.5|) which then yields the solution (jl.2j) for the recurrence (|1.1|) . It is convenient for us 
to summarize the results of [21] by expressing the solution of this initial value problem in 
the form of an algorithm, as follows: 

Step 1: Find the backwards iterate r_i from the initial data, evaluate the quantities /o — 
r 1 r_ 1 /(r ) 2 and f\ = t 2 t /(ti) 2 , and use these to calculate the integral 

J = J(foJi) = p"(k) 



Step 2: Use J, a, (3 to calculate 

This gives the point P = (A,/i) = (p(n), p'(k)) e E, with /i = ±^/a. 

Step 3: Construct the invariants g 2 , <?3 of the curve E as in (jl.4|) . from the formulae 

g 2 = 12A 2 - 2J, g 3 = 4A 3 - g 2 X - a. 



Step 4: Iterate (jl.5|) backwards to obtain /_i from and j\. Hence find the point Pq = 
(u, £) = (p(zo), p'(zq)) e E from the formulae 



\ f c /o CA ~ f-i) 

v — A — Jo, ? — • 

Step 5: Calculate the values k, z ^ C from the elliptic integrals 

/• p dx f p ° dx 

k = ± — , z = ± — ; 

Joe y Joo y 

these should be interpreted as the points in the Jacobian Jac(E) corresponding to the points 
P,Po £ E respectively. Note that because of the involution y — > — y these values are only 
defined by to an overall ± sign, subject to the constraint that p'(k)p'(zq) = £/i = /f (fi~ f-i) 
as in Step 4. Once zq and k are obtained then A are B are found from the formulae (jl.3|) . 
Remarks. It is useful to note that the coefficients a,/3 in the recurrence are given as 
elliptic functions of k by 



a 



p'(K) 2 , (3 = p\k) 2 {p{2k)-p{k)Y (2.8) 



The above solution of the initial value problem establishes an exact correspondence between 
two sets of six parameters: the parameters gi, g%, k, Zq, A, B that specify the elliptic curve E, 
the two points P, P G E, and the prefactors in (|1.2j) : and the parameters a, (3, r , Ti, t 2 , t 3 
specifying the constant coefficients and initial data for the recurrence (J1.1J) . In order to 
interpret (jl.5)l as an integrable map, it is necessary to further specify a symplectic structure 
0113]; symplectic coordinates and a Lax pair were given in j23], which make equivalent 
to the discrete g — 1 odd Mumford system in |28j . 

As an example of the above algorithm, we present the results for the Somos (4) sequence 
(|2.7|) . with a = (3 = r = T\ = r 2 = r 3 = 1. We find r_i = 2, so /o = 2, /i = 1 gives J = 4 in 
Step 1. In Steps 2 and 3 we have A = 1, set fi — 1 and then find #2 = 4, #3 = — 1, and in Step 
4 we obtain /_i = 3/4 so that ^ = — 1, £ = 1. Thus the Somos (4) sequence corresponds to 
the sequence of points P + nP = (—1, 1) + n(l, 1) on the curve 

E : y 2 = 4x 3 - 4x + 1. 

Finally, evaluating the elliptic integrals and sigma functions to 9 decimal places using the 
MAPLE computer algebra package (version 8), we find that the curve has real and imaginary 
half-periods uj\ = 1.496729323 and = 1.225694691z respectively, while 

/oo 
(At 3 - At + dt = -1.134273216, 

/oo 
(4t 3 - At + dt = 0.204680500 - 1.225694691z, 

which yield the other quantities in ()1.2j) as 

<t(k) = 1.555836426, A = — = 0.112724016 - 0.824911687i 

B = ) ' y °( = 0.215971963 + 0.6160281931 
a(z + k) 



However, the sequence of arguments of the sigma function can be written more succinctly as 

z + nn = (2n - 3)£b, 5 = 0.929592715 + u 3 , 

so that the iterates of the recurrence correspond to the sequence of points (2n — 3)P on the 
curve E, where P = (0, 1), P = 2P. The full sequence of points nP is associated with the 
elliptic divisibility sequence (12. lj) . 

Elliptic divisibility sequences are currently of considerable interest due to the fact that 
large prime numbers can occur therein (i.e. t p may be prime when the index p is prime, see 
O [13 E]). Cantor has considered the division polynomials for odd hyperelliptic curves [TT] . 
corresponding to sequences of divisors n(P — oo), which also satisfy higher order recurrences 
written in terms of Hankel determinants; Matsutani has obtained the functional form of 
these division polynomials in genus two jHO]- In the next section we shall derive an eighth 
order bilinear recurrence associated with the sequence of divisors T> Q + n(P — oo), where 
T>q = (Pi — oo) + (P 2 — oo) is the reduced divisor of two points on a hyperelliptic curve of 
genus two. 



3 Addition of one point in genus two 

Let us consider an algebraic curve X of genus two defined by the affine model 



X :-- 



v 



4 

f(x) =4x 5 + ^ Ci x J }, 



(3.1; 



j=0 



which realizes the curve as a two-sheeted covering of the Riemann sphere with 2g+l branch 
points in the complex plane plus a single branch point oo at infinity. The vectors of canonical 
holomorphic differentials and canonical meromorphic (second kind) differentials are denoted 



du 













dr = 


If ) 









/ (12x 3 +2c 4 x 2 +c 3 x)dx \ 
4y 



V 



x z dx 

y 



J 



respectively. If we let (A±, A 2 ; B\, B 2 ) denote the canonical homology basis for the compact 
Riemann surface corresponding to X, with non-vanishing intersections Aj ■ = 5jk, then 
the 2x2 matrices of A- and 5-periods are given by 



9a 



2uj 



L - 

J Ax y JA2 y 

£ x dx f x dx 

JAi y JAi y 



£ dx 

JE 



r dx 

JBx y JB2 y 

£ x dx £ x dx 

JBi y JB2 y 



The Jacobian of X is the complex torus Jac(X) = C 2 /r, where T = 2u7? © 2uj"I? is 
the lattice generated by the periods of canonical holomorphic differentials. The elements 
(Pi,P 2 ) of the symmetric product (X) 2 can be identified with degree zero divisors T> = 
(Pi — 00) + (P 2 — 00), which are mapped to Jac(X) by the Abel map: 



u 



P-2 



du + / du 6 Jacpf ) 



(where here we are basing the map at 00). 



The Kleinian sigma function cr(u), which is a quasiperiodic function of (u\, w 2 ) T = u 6 C 2 , 
is the genus two analogue of the Weierstrass sigma function. The Kleinian £ and p functions 
are defined by 

dlogo-(u) . 
0(u) = ^— , J = 1,2, 

. . <9 2 loga-(u) . . (9 3 logo-(u) . , „ n 

We refer the reader to other works such as [3 IE] for a detailed introduction to hyperelliptic 
curves, Kleinian functions and their definition in terms of Riemann theta functions (see also 
jH H21 ESI El and references). 

Vectors v in the theta divisor in Jac(X) can be characterized by the fact that 

o-(v) = 0. 

We wish to take a vector v 6 6 in the theta divisor, given by 

du; (3.2) 

so v corresponds to the (reduced) divisor of a single point (A, /x) G X, with /x 2 = /(A). 
In |Sj (Theorem 4.9) it is proved that the Baker function 

$ : Jac(X) xI^C, 

defined for u e Jac(X) and (A, /x) G X by 



$(u; (A,/x)) = J V x /^n ex P I -f dr r u ) , v=/ du, 

satisfies the Schrodinger equation 

(<9 2 -2p 22 )$ = (A + c 4 /4)$. (3.4) 

Note that we have chosen a particular normalization for the Baker function compared with 
jH], including the denominator 02 (v) 7^ 0, and the principal value symbol § in ()3.3|) denotes 
the fact that the integral of the meromorphic differential dr is regularized at infinity. 

Let us define two different Baker-Akhiezer functions $>± related by the hyperelliptic 
involution, as 

$± = $(u; (A,±/x)). 
Then from the proof of Theorem 4.9 in jS] (restricting to g — 2) we have that 

o 1 ^ ±/x + <9 2 P(A; u) 

a 2 log<l> :t = — — — ^ — -, (3.5) 

& 2P(A;u) ' v ; 

where 

V(X; u) = A 2 - p 22 (u)A - p 12 (u) (3.6) 

is the Bolza polynomial in genus two [S]. Hence it follows that V = V(X;u) satisfies the 
Ermakov-Pinney equation with respect to derivatives <9 2 in the variable xt 2 , namely 

1 u 2 
V(d 2 2 V)--(d 2 V) 2 + 2VV 2 + ^-, V := -2p 22 (u)- A-C4/4. (3.7) 



It is a well known classical result of Ermakov (see [20] for references) that the general 
solution of the Ermakov-Pinney equation ()3.7j) is just given by a product of two solutions of 
the Schrodinger equation (d 2 + V)ip = with Wronskian ^. 
Taking the difference of the ± equations (|3.5|) we have 

_ /i 

where 

cf> $ . 
)= ~ + 
v +/ <9 2 $_ <9 2 $+ 

is the Wronskian. Clearly this must be independent of u 2 , but we claim that in fact this 
Wronskian has precisely the value /i, which means that the product $ + $_ = V. Rewriting 
this in terms of the sigma function, we can state the following result. 

Proposition. The Kleinian sigma function for a hyperelliptic curve \3. 1)) of genus two 
satisfies the following formula for addition of a single point on the curve: 

.(u + vMu-v) 

cr(u) 2 o- 2 (v) 2 

In the above, u G Jac(X) is a generic vector in the Jacobian, vG0C Jac(X) is the image 
of the single point (A, /i) G X under the Abel map, and V is the Bolza polynomial defined by 
(EU). 

Proof: Starting from Baker's addition theorem for genus two [2], 
cr(u + v)er(u — v) 

- ( v2 , , 2 = P22 (u )pi2 v - p 12 u )p 22 v + p n v - pii(u), 3.9 

where u, v are generic points in Jac(X), and multiplying both sides by a(\) 2 /a 2 ('v) 2 , the 
result follows by taking the limit cr(v) — > as v tends to the theta divisor. It is necessary 
to use the fact (see e.g. [22]) that the x coordinate of the point (A, /S) 6 X is given, in terms 
of derivatives of the sigma function evaluated on the theta divisor 0, by the expression 

A = -^4 



This follows from the fact that differentiating er(v) = with respect to A gives, by the chain 
rule, 

dvi dv 2 o-i (v) Acr 2 (v) 

— dia v + — d 2 a v = + = 0, 

d\ aX /i /i 

for v G 6 given by (|3.2j) . □ 

Remark. Enolskii and Gibbons recently calculated the exact analogue of (|3.8|) in genus 
three. The addition formula (J3.8|) is a special case of the generalized Frobenius-Stickelberger 
addition formula in genus two considered in 13*3*] . Onishi has further generalized the 
Frobenius-Stickelberger formula to hyperelliptic sigma functions for all genera [21], and the 
special case of the formula corresponding to addition of one point has been applied to the 
problem of construction of Wannier functions for quasi-periodic finite-gap potentials in [1] . 

Cantor has constructed the division polynomials for hyperelliptic curves, and obtained 
certain recurrence relations for them in the paper where in particular an eighth order 
bilinear recurrence is found in genus two. Up to a suitable normalization, Matsutani has 



considered the exact analytic expression for these division polynomials, which are equivalent 
to the sequence of functions 

a(nv) 



a 



known as as hyperelliptic psi-functions flty I3()|. In the following theorem, we present a se- 
quence of tau-functions that generalize these psi-functions and yet satisfy the same recurrence 
of Somos 8 type. 

Theorem 2. Define the sequence {r n \ n G Z} by 

r n = AB"^±^l, (3.10) 

cr 2 (v) n 

where u G Jac(X) is a generic vector in the Jacobian of the genus two curve v G 

O C Jac(X) is the image of the single point (A,/i) G X under the Abel map, a denotes the 
Kleinian sigma function of the curve, and A, B are arbitrary constants. Then the terms of 
the sequence satisfy a bilinear recurrence of order 8, given by 

3 

T~n+A T n-& = ^ ] Ctj r n+j r n-j, (3- 11) 

j=0 

where the coefficients ctj ( independent of n) are given by 

^(6vM3v) 2 a(3v)a(5v) 
1 a(4v)a(2v)V 2 (v) 30 ' 3 a(2v)a(4v) ( x 2 (v) 14 ' 1 ' ' 



q(4v) 2 / _ q(3v) 3 q(5v) 



cr 2 (v) 24 <T(2v)^ \ <j(4v) 3 a(2v) J 

Proof. Substituting the expression (j3.10j) into (j3.11|) and using Baker's formula (|3.9J) 
together with the result of the Proposition yields an expression of the form 

C (v) + Cn(v)pn(u) + Ci 2 (v)pi 2 (u) + C 22 (v)p 22 (u) = 0. 

The three functions pjk(u), j, k = 1, 2 on Jac(X) are not linearly dependent (although they 
do satisfy a nonlinear relation [HIE]; giving the Kummer surface in CP 3 ). Therefore each 
of the coefficients Co(v), Cn(v), C*i 2 (v), C 22 (v) must vanish, which leads to a linear system 
for the ctj as functions of v. This determines the above formulae for the coefficients ctj 
uniquely in terms of a, its derivative <r 2 , and the pjk evaluated at various multiples of v. 
The terms involving can be removed by making use of the addition formulae (J3.9)) and 
(J3.8|) to yield the expressions (|3.12j) and (J3.13)) in terms of a and <r 2 alone. Taking the limit 
u — > 0, these are equivalent to Matsutani's expressions for the coefficients in the eighth order 
bilinear recurrrence for the psi-function (see formula (3.13) in j2H])- 1=1 
Corollary. The sequence of Bolza polynomials 

fn = V{\\ u + nv) = A 2 - p 22 (u + nv)A - pi 2 (u + nv), (3.14) 

for u G Jac(X) and v G 0, satisfies the sixth order nonlinear difference equation 

3 3 j-l 

fnWiUkfn-k)^ = «° + Y. a 'jfnX\Un + kfn- k y-\ (3.15) 
k=l j=l k=l 

with the coefficients ctj as given in equations and \3.1§) . 



Proof of Corollary. Upon setting 

with r n as in ()3.10|) . and using the addition formula ()3.8j) . the result is an immediate conse- 
quence of Theorem 2. □ 

Remarks. Taking u as the Abelian image of (Pi,P2) G {X) 2 , the sequence f!3.14|) corre- 
sponds to the linear flow u + nv in the Jacobian, or equivalently the sequence of divisors 
V n ~ £> + n(P - oo) ~ V - n(P - oo) with V ~ (P 1 - oo) + (P 2 - oo) and P = (A, -pi), 
the image of P = (A, fi) under the hyperelliptic involution. The Bolza polynomial leads to 
the solution of the Jacobi inversion problem for the curve (J3.1|) (see Theorem 2.2 in jHj, and 
section 4 below), so that in particular if V n ~ (xi(n), yi(n)) + (a; 2 (?7,), ?/ 2 (n)) — 2oo, then we 
have 



f n = P(A; u + nv) = (A - xi(n))(A - x 2 (n)), yj(n) = -<9 2 P(A; u + nv) 



\\=Xi 



1,2. 



Cantor's results in concern the sequence of reduced divisors T> n ~ n(P — oo), correspond- 
ing to the multiples of a single point on an odd hyperelliptic curve of genus g. In particular 
for g — 2 he obtains a bilinear recurrence of order 8, which is the degenerate case u = 
(Pi — > oo, P 2 — > oo) of our construction, while the analytic derivation in that case appears 
in the work of Matsutani ED]- The sixth order difference equation (|3.15|) appears as 
equation (3.15) in [2H], where the special solutions with u = are also presented. 

The sigma functions of genus g odd hyperelliptic curves, given by y 2 = f(x) with / a 
polynomial of odd degree 2g + 1, are known to be tau functions of the Korteweg-deVries 
(KdV) hierarchy of partial differential equations (see jS] for instance). It is also known that 
when the curve degenerates completely to y 2 = 4x 2g+1 , the corresponding sigma function 
degenerates to a polynomial (see jUIHl]); which gives a rational solution of KdV in terms of 
a Schur function (see £Q and chapter 14 in It is instructive to consider the case when 

the curve (|3.1|) for g = 2 degenerates to a singular rational curve: 



y 



Ax 5 . 



In that case the Kleinian sigma function degenerates to the Schur function 



o u 



U 2 
Ul , 

3 ' 



(3.16) 



which is the tau function of the three-pole rational solution of the KdV equation 

49iU = dlV + 6Vd 2 V, V = -2p 22 (u) = 2d\ logtr(u). 
The theta divisor consists of vectors of the form 



/ K-M Ax \ 



u 



4a: 5 
x Ax 



lA- 3 /2 

3 



-A-V2 



7 3 /3 
7 



7eC, 



satisfying cr(v) = 0. It is trivial to check that the Schur function satisfies the addition 
formula (J3.8|) . Defining r n in terms of the Schur function (J3.16j) by (J3.1U)) . it is easy to verify 
that this gives a particular solution of the eighth order recurrence (|3.11|) with 
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4 BT for the case (ii) Henon-Heiles system 

The integrable case (ii) Henon-Heiles system is a system of two degrees of freedom defined 
by the natural Hamiltonian 



2 



777/ 

hi = ^ (Pi + Pi) + 4 + 2«2?i _ 2 aq ^ + Cq2 ~ 2gf ■ ^'^ 
The coordinates qj and momenta pj are canonically conjugate, and Hamilton's equations 

^ = {h u q 3 h ^ = {h u p,}, j = 1,2 (4.2) 

are equivalent to the ordinary differential equation for travelling wave solutions of the fifth 
order flow in the KdV hierarchy [TJj. The equations of motion ()4.2j) can be written in the 
form of a Lax equation 

dL 



where the Lax matrix L is 



L(x; qj , Pj ) = ( _ 2 jf X i 2 f ( ^ S'^ ) > ( 4 - 3 ) 

2ar + (g 2 + 2a)x - ±q( -±p 2 x + |pi<?i 1 



with 



B(x] qj,pj) = 2x 3 + (-q 2 + 4a)x 2 + 1 (qi + Ml - aq 2 + 16a 2 + 4c) x + \ (p\ - 7 ^ T ] . 

o 8 V <?i / 

(Note that we have L — > — '■""^ L, it — > x + a compared with reference |2*T|.) The Lax 
equation is the compatibility condition for the linear system 

L\I/ = yty, — = N\l>, N = ( j o J ■ ( ^ 

The genus two spectral curve is of the precise form (|3.1|) . namely 

777, 

det (L - yl) = y 2 - 4a; 5 - 12ax 4 - (c + 12a 2 )x 3 - -M 2 - -h 2 x - — = 0, 



with 



1 2 1 1,1, 2 ,2\2 rn?{q 2 + 2a) 

^2 = 7W2 + 2a)p! - -gipip 2 - - o (?2 - 2a <?2 + c + 4a ) ^ 



2 



being the second independent integral, in involution with fai i.e. {hi,h 2 } = 0. The integral 
h 2 generates a second commuting flow 

^ = {h 2 , g] }, ^ = {h 2 , Pj }, j = 1,2. 

Up to a shift of origin, the time variables s, t can be identified with the coordinates 
U\,u 2 respectively on Jac(X). Using the results of Theorem 2.2 in jH], the solution of the 
Henon-Heiles system can be reduced to the Jacobi inversion problem 



u = / du + / du, 



where the separation coordinates x±,x 2 are found from the (2,1) entry in the Lax matrix 
(|4.H|) . given as a multiple of the Bolza polynomial by 

2x 2 + (q 2 + 2a)x - \q\ = 2{x - x x ){x - x 2 ) = 2V(x; u). 
o 

The separation variables (xj,yj), j = 1,2 correspond to the reduced divisor V = (xi,yi) + 
(x 2 ,y 2 ) — 2oo, and they are related to the Kleinian functions by 

Pl2(Uj = -XxX 2 , p22(u) = Xi + X 2 , p22l(Uj = , p222(» = 

X\ — X 2 Xi — x 2 

(see e.g. [12] )• Thus the connection with the Bolza polynomial immediately leads to the 
solution of the Henon-Heiles system in terms of Kleinian p functions, which is 

q\ = 16pi 2 (u), q 2 = -2p 22 (u) - 2a, qipx = 8p 22 i(u), p 2 = -2p 2 22(u). 

It is shown in |2I] that the case (ii) Henon-Heiles system has a Backhand transformation 
(BT) with parameter A, which is a symplectic map with generating function F(q\, q\\ A) such 
that dF = J2j=i 2 Pjdqj — Pjdqj- The explicit form of the generating function is 

F = Z + | log (|^) + f^ 5 + 4(g 2 + q2)Y 3 + (2(q 2 + ~q 2 f - 2q 2 q 2 + \{ql + q\) + 2c) F, 
where Z(qj,qj) and Y(qj,qj) are defined by 

Z 2 = m 2 + Agfo 2 , F 2 = A + a - - (q 2 + q 2 ) . 
The BT can be realized as a similarity transformation on the Lax pair (discrete Lax equation) 

LM = LM, (4.5) 

where L = L(x; qj,pj) and 

is the elementary Darboux matrix (see [HI]). Clearly from ()4.5|) the BT preserves the spec- 
trum of L, and so maps solutions to solutions. 

In fact, the BT was constructed in [21] by making use of the formulae for the Darboux 
transformation of the Schrodinger equation, since the components of ^> in the linear system 
()4.4|) are given by \I/ = (i/j, d 2 ^) T with (d^ + q^ip — (x + a)ip. Then the quantity Y appearing 
in the Darboux matrix can be given explicitly in terms of the Baker function $ + defined in 
(jH3D as 

Y = d 2 log $+ = d 2 log $(u; A) = C2(u - v) - £ 2 (u) + -f dr 2 , 

J oo 

and by a simple calculation using ()3.5|) we also have 

ti + d 2 V(X; u) -n + d 2 V(X; u - v) 



2P(A;u) 2P(A;u-v) 



F 

which is equivalent to equation (2.9) in 

It follows from general results in [21] that the BT defined by the Darboux matrix ()4.6|) 
just gives a shift on the Jacobian, and here we see that 

qf = 16pi 2 (u- v), q 2 = -2p 22 (u- v) -2a, gipi = 8p 22 i(u - v), p 2 = -2p 222 (u - v). 



In terms of divisors, we have the equivalence 



V ~ V + (A, -/i) - oo ~ (xi, yi) + (x 2 , y 2 ) - 2oo ~ V-\. 

Similarly, applying the same BT but using the Baker function $_ corresponds to adding the 
point P = (A,/i) instead of P — (A, —fi), which gives the divisor T>\ ~ V + (A,//) — oo. 
In other words, the nth term in the sequence of Bolza polynomials (j3.14j) . with argument 
u + nv G Jac(X), is the result of n applications of the Henon-Heiles BT, adding the same 
point P each time. We should point out that the Henon-Heiles BT given here and in [21] 
is the same as the g = 2 discrete odd Mumford system in [2H], with a particular choice of 
symplectic structure. 

5 Conclusions 

We have constructed an eighth order bilinear recurrence relation (|3.11|) whose nth term is 
expressed in terms of the Kleinian sigma function of an odd hyperelliptic curve (|3.1|) of genus 
two. Moreover, we have explained how this is connected to the BT for the case (ii) integrable 
Henon-Heiles system of |2T| 122] : each shift n — > n + 1 in the recurrence corresponds to the 
shift on the Jacobian induced by the BT. 

The question remains of whether we can solve the initial value problem for an eighth 
order recurrence of the form (|3.11|) . In fact, the expression ()3.10|) depends on at most 10 
parameters: the five coefficients Cj, j = 0, . . . , 4 that specify the curve plus the three points 
Pi, P 2 , P on the curve to specify the sequence of divisors V n ~ (P 1 — oo) + (P2—oo)+n(P—oo), 
and the two prefactors A, B. (Actually, the constant C4 can be removed from the start by 
making a shift in x.) The general eighth order recurrence (J3.ll)) would have four parameters 
otj, j = 0, 1, 2, 3 in general position and 8 initial data Tj, j — 0, . . . , 7. In fact, by counting 
arguments, we expect that the general solution of a bilinear recurrence consisting of N + 2 
bilinear terms, 

N 

T n+N+l T n-N~l = ^ a j T n+j T n-j, 
3=0 

which is of order 2N + 2, should correspond to a sequence of divisors V n ~ ^2f=i{Pj ~ 00 ) + 
n(P — 00) on an hyperelliptic curve of genus N, since this solution should depend on 3N + 3 
parameters. This expectation agrees with the form of vector addition formulae for Riemann 
theta functions obtained by Buchstaber and Krichever [7j, which have precisely iV-|-2 terms 
in genus N. However, preliminary calculations suggest the need to consider a different model 
for these curves, without a branch point at infinity. For example, sigma functions for curves 
of the form 

5 

y 2 = 4x 6 + ^2 C J xj 
3=0 

should give solutions to the Somos 6 recurrence 

2 

T~n+3^~n— 3 ^ ^ OCj T n +jT n —jj 
3=0 

this should correspond to the BT for the g = 2 even Mumford system in [2H], which has a 
spectral curve of this type. This agrees with recent results of van der Poorten which 
show that a certain class of Somos 6 sequences arise from the continued fraction expansion 
of the square root of a sextic. 



Note added in proof. In a private communication, David Cantor has shown one of us 
that by extending the results of [TT] to divisor sequences V + n(P — oo) it is possible to 
prove that the sequence of tau functions r n given by ([3.10]) in Theorem 2 satisfies a family 
of trilinear recurrence relations given in terms of a Hankel type determinant, that is 



m ' n 1 n+mTn—m 



0"m T n-2 
0>m—l7~n—l 
Q"m—2l~n 



Om-lTn+1 O'm'Tn+2 



(5.1; 



for all m, n G Z, where 



(71 mv 



0~2 (V 



The above formula is the genus two analogue of a family of higher recurrences satisfied 
by Somos 4 sequences j^HI (cf. also Morgan Ward's identity (|2.6j) for elliptic divisibility 
sequences). Up to a suitable choice of normalization, these a m are the genus two division 
polynomials derived by Cantor [TJ, whose functional form has been specified precisely in 
the works |2HI3U] of Matsutani, where they are referred to as hyperellliptic psi-functions. In 
the limit u — > 0, when r n — > a n , the equation ()5.1|) reduces to a Hankel formula for these 
psi-functions obtained by algebraic means in [11] . which is also presented as formula (3.10) 
in and rederived using analytic techniques in |3U]. Setting m = 3 in ([5.1]) yields a 
sixth order trilinear recurrence for r n , while setting m = 4 yields an eighth order trilinear 
recurrence. As was first pointed out to us by Cantor, the bilinear recurrence ([3.11]) above 
then follows as a consequence of these two trilinear identities, upon eliminating between the 
two of them. Interestingly, Buchstaber and Leykin have very recently derived a trilinear 
differential addition formula for genus two sigma functions - see Theorem 5.6 in [TU] . 
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